clear;
T=10;
N=10000;
dt=T/N;%diataT
t=0:dt:T;%时间
y0=8;%初始条件
y1=zeros(1,length(t));%建立一个全是0的向量
y1(1)=y0;%将初始条件放入
v0=8;%初始条件
y2=zeros(1,length(t));%建立一个全是0的向量
y2(1)=v0;%将初始条件放入
for n=1:length(t)-1
    y1(n+1)=y1(n)+y1(n)*dt;
    %y2(n+1)=y2(n)+(-k/m*y1(n))*dt;
    y2(n+1)=y2(n)+y2(n)*dt+1/2*y1(n)*dt^2;
end
%plot(t,y1,t,y2);
plot(t,y2,t,y1)
%利用holdon更改划分，判断怎么样的精度合适
hold on